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(~| . Abstract 
^ ■ 

We provide (partial) reconstruction formulas and discrete Fourier transforms for wave 
functions in standard Fock-Bargmann (holomorphic) phase-number representation from 
a finite number N of phase samples {9k = 2'n:k / N}^^^ for a given mean number p of 
\ particles. The resulting Coherent State (CS) subsystem S — {\zk = p^e^^'')} is complete 

^ I ■ (a frame) for truncated Hilbert spaces (finite number of particles) and reconstruction 

I formulas are exact. For an unbounded number of particles, S is "almost complete" (a 

pseudo-frame) and partial reconstruction formulas are provided along with an study of 
OA I the accuracy of the approximation, which tends to be exact when p < N and/or N ^ oo. 

O : MSG: 81R30, 81R05, 42B05, 30H20, 42C15 

Keywords: coherent states, discrete frames, reconstruction formulas. 



1 



1 Introduction 



The subject of completeness of Coherent State (CS) subsystems traces back to von Neumann 
results concerning the density of canonical CS [Ij and is closely related to standard theorems like 
Fourier, Heisenberg (uncertainty relations), Nyquist-Shanon (information theory), etc. Criteria 
of completeness for a CS subsystem S = {\zk)} corresponding to a set of points Q = {zfc} 
in the complex plane C are known in the literature (see e.g. the standard text books ^ i3j 
and references therein). In general, S is complete iff no function (vector symbol) ^{z) = {z\^l'), 
7^ 0, in the Fock-Bargmann space, vanishing at all points in Q, exists. Of special interest is the 
case where the points z^ form a rectangular lattice z^ = Zmn = rnui + nu}2 , where the periods 
of the lattice wi and UJ2 are linearly independent, S = 9(a;2Wi) 7^ 0, and m,n are arbitrary 
integers. The simplest example is the von Neumann square lattice with elementary Planck cell 
area S = ir, corresponding to the most accurate simultaneous measurement of both coordinate 
and momentum. For S = n, the CS subsystem S = {\zmn)} is complete and it remains complete 
if one state is removed, but becomes incomplete when any two states are removed. For S < n, 
the CS subsystem S is overcomplete and remains so if a finite number of states are removed. For 
S > TT, the CS subsystem S is not complete. This gives a characterization of overcompleteness 
in terms of density of states per Planck cell. All these results can also be stated in the language 
of (discrete) frames and sampling and interpolating sets (see for instance [4]), and we shall make 
use of this terminology in this article. 

The square lattice Q = {zmn} provides a sampling of the (position-momentum) phase space, 
so that reconstruction formulas for arbitrary functions ^{z) can be obtained from its samples 
^(zmn) for S < TT. In this article we shall use an alternative phase-number parametrization of 
z = re*^ in terms of "mean number of particles" p = r"^ = |zp and phase 6 = arctan(3'(z)/3?(2:)) 
and we shall choose a CS subsystem S = {\zk)} corresponding to a finite set of N points Q = 
{zk = re^'^'''^f^}^~Q uniformly distributed on a circle of fixed radius r. This would physically 
correspond to a finite number N of phases {Ok = ^'^'''k/N}^~Q for a fixed mean number of 
particles p = r"^. Of course, the resulting CS subsystem S = {\zk)} won't be complete for finite 
N but partial (almost exact) reconstructions of functions ^{z) from samples {^{z^)} will 
be possible and we shall study the accuracy of the approximation. We shall also study the 
conditions under which an exact reconstruction of ip is possible, namely when restricting oneself 
to a truncated Hilbert subspace (finite number of particles) or for mean number of particles p 
below the critical value pc = N. The Fourier coefficients (or probability of occupation numbers) 
can be obtained by means of the (filtered) Discrete Fourier transform of the data, allowing for 
a straightforward fast extension of the reconstruction algorithm. 

The organization of the paper is as follows. In Section [2] we briefiy review some properties 
of canonical CS and their restriction to the circle, suggesting the possibility of reconstructing 
states from phase samples for a finite mean number of particles. Section [3] discusses the case of 
exact reconstruction for truncated (finite-dimensional) Hilbert spaces (finite number of particles- 
bosons) in the language of discrete frames, sampling operators and "sinc-type" kernels. The 
more interesting and general case of partial reconstruction of arbitrary functions in terms of 
finite (non-complete) CS subsystems (pseudo- frames and undersampling), in the number-phase 
parametrization, is developed in Section |H The advantage of this number-phase parametrization 
is that it allows the explicit inversion of resolution and reproducing kernel operators by making 
use of the theory of Circulant Matrices and Rectangular Fourier Matrices. We also study the 
accuracy of the approximation and the conditions under which it tends to be exact. Section [5] is 
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left for some comments on potential physical applications in Quantum Optics and an outlook. 
We also introduce a couple of appendices: Appendix A commenting on traditional problems 
associated with the proper definition of a phase operator, and Appendix B with the proof of the 
main theorem of our work. 



2 Canonical Coherent States and its Restriction to the Circle 

Bose quanta are described by a pair of annihilation (lowering a) and creation (rising a^) ladder 
operators fulfilling the commutation relations: 

[a,a^] = l. (1) 

The Hilbert (occupation-number or Fock) space H realizing this commutation relation is made 
of (normalized) eigenstates |n) of the (Hermitian) number operator J\f = a^a, i.e.: 

M\n) = n\n), {n\m) = 5nm- (2) 

These states can be generated from the Fock vacuum |0) as: 

\n) = -^(atno) (3) 



and they fulfill the completeness condition X^J^o ~ ^■ 

Coherent states z £ C, are defined as eigenstates of the annihilation operator a, 

a\z) = z\z). (4) 

They can be explicitly generated by acting with the (unitary) displacement operator 

U{z, z) = e^"^-'"'' = e-"""/2e^'^^e^^ (5) 

on the vacuum |0), that is: 

oo ^ 

\z) = U{z,z)\{)) = e-^^/2gW|Q^ =e-^^"/2 (6) 

Since a is not hermitian, there is no reason (in principle) for the set {Iz),^; G C} to span 
an orthonormal complete set. In fact, one can easily compute the CS overlap (or Reproducing 
Kernel) 

C{z, z') = {z\z') = e-^^"/2e-^'2'/2g5^'^ (7) 

which is not zero for z 7^ z' . Actually, the set {\z), 2; G C} turns out to be overcomplete. In fact, 
it defines a tight frame in % , with resolution of unity 

I= - f \z){z\(fz, (8) 



vr _ ^ 

where we denote d?z = (iRe(z)c?Im(z) = rdrdO the Lebesgue measure in C in polar coordinates 
z = re*^. Indeed, using the expansion ([6|) we have that 



- / \z){z\(fz = - / e-^'^ y 4^ 
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-1 ^ />oo „n+m+i i-z-K 

- y / dr , |n)(m| / d^e^^^""™) 

n,m=0 

= 2 y dr j — e"*" |n)(n| = |n)(n| = /. 



n=0 n.=0 

The average number of particles p (bosons) in a coherent state \z) is 

p = {z\Ar\z) = \z\^ (9) 
and the probabihty ampUtude of finding n particles in \z) is: 

Uniz) = {n\z) = {n\U{z,z)\0) = e-'~'/^-^z'^. (10) 

vn! 

In the same way, the probability amplitude of finding any state 

oo 

\i;) = Y,(^n\n)en (11) 

n=0 

in the coherent state \z) is given by: 

oo 

^{z) = {z\i;) = Y,aJUz). (12) 

n=0 

The "Fourier" coefficients can be calculated through the following integral formula: 
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an = (n|V^) = - / ^{z)Un{z)d'z. (13) 

Due to the over-completeness of the coherent states and to the holomorphic character of the 
representation, the integration on the radius |z| is not really necessary. Actually, we can recover 
a-n from the values of ^{z) on z = ^/pe*^ for a fixed average number of particles p = |zp as: 



/ n! 1 

p^e-P 27r 



ei"«^r(^e*'^)d0. (14) 



This fact suggests the possibility of reconstructing the state ip from a finite number A'" of phases 
{^k}k=o ^^'^ ^ fixed (but arbitrary) mean number of particles p. We shall pursue this idea in the 
next Sections and we shall study the conditions under which total and partial reconstructions of 
ip from N samples {^'(zfc), = \/P^*^'°}£o^ possible and the accuracy of the approximation 
for partial reconstructions, which tends to be exact in the limit N — )• oo. 

We shall not deal here with the traditional interesting problems associated with the proper 
definition of phase operators (see e.g. O [U [71 [51 [9l [lOl E]). Instead, we address the reader 
to Appendix [X] for a short review on the subject. We shall just point out that the inherent 
difficulties in defining a proper phase operator can be overcome by restricting oneself to a 
truncated Hilbert subspace T-Lm of T-L (finite number of particles). One can also circumvent 
this problem by considering quantum systems with a large average number of particles (i.e., 
semi-classical states), for which phase (and number) uncertainties are negligible. 
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3 Finite Number of Particles 



Firstly, we shall work with systems with finite (but arbitrarily large) number of particles M < oo, 
so that our Hilbert space will be the truncated state space Hm T~i of M + 1 dimensions. 
Following the terminology introduced in [HI |9] , we denote by physically accessible o preparable 
states as those that can be obtained from the vacuum state by acting during a finite time with 
a source of finite energy and with a finite interaction. These states belong to Hm for some 
integer M. There are also states, denoted by physical states, that are characterized by the fact 
that {M'^) is finite for finite, but arbitrary q £ N. These states include coherent states, squeezed 
states, thermal equilibrium states and in general Gaussian states [12]. 

Before discussing the reconstruction of states ip £ T-Lm from a finite number of phase samples, 
we shall remind some basic definitions on discrete frames, sampling operators and sinc-type 
kernels. The advantages of using finite tight frames for finite-dimensional Hilbert spaces has 
been recently discussed in [T^. 

3.1 General considerations on discrete frames 

Our purpose is to discretize the integral ([8]), by restricting ourselves to a discrete subset Q = 
{zk} C C. The question is whether this restriction will imply a loss of information, that is, 
whether the set S = {\zk),Zk G Q} constitutes a discrete frame itself, with resolution operator 

^ = Yl \^k){zk\- (15) 

^feGQ 

The operator A won't be in general a multiple of the identity. In fact, a continuous tight frame 
might contain discrete non-tight frames, as happens in our case (see later on Sec. I3.2p . 

Let us assume that S generates a discrete frame, that is, there are two positive constants 
< & < -B < oo (frame bounds) such that the admissibility condition 

^^iiv'f < E \{^km\'<Buf (16) 

holds \/tp G T-Lm- To discuss the properties of a frame, it is convenient to define the frame 
(or sampling) operator T : H-m ^ ^"^ given by T('0) = {{zk\'4^), Zk G Q}- Then we can write 
A = T*T, and the admissibility condition (jl6p now adopts the form 

bl < T*T < BI, (17) 

where / denotes the identity operator in "Hm- This implies that A is invertible. If we define the 
dual frame {\z) = A^^\z)}, one can easily prove that the expansion (reconstruction formula) 

IV') = E ^k\~Zk), (18) 

where = (-^fclV')) holds in Hm, that is, the expression 

T+T = \zk){zk\ = T*{T+r = Y \^k){zk\ = I (19) 

Zk^Q ZkGQ 
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provides a resolution of the identity, where 7^ = (7'*7')-i7'* ig the (left) pseudoinverse (see, 
for instance, |14] ) of T (see e.g. [3l [15] for a proof, where the authors introduce the dual frame 
operator T = (7^^)* instead). 

It is interesting to note that the operator P = TTf^ acting on is an orthogonal projector 
onto the range of T- Eventually, we are interested in a finite number N of phase samples, so 
that the space i"^ can be substituted by C^, with N the cardinal of Q. 

Prom (lisp the function ^{z) can be obtained 

^{z) = {z\^l^) = Ek{z)^k (20) 

from its samples = (zkli^), through some "sinc-type" kernel 

Ek{z) = {z\zk) (21) 

fulfilling Hfc(z/) = Pik- A projector is obtained, instead of the identity, to account for the fact 
that an arbitrary set of overcomplete data G C^, can be incompatible with E Hm, 
and therefore they are previously projected (note that an overdetermined system of equations 
is being solved). 

This case will be named over sampling, since there are more data than unknowns, and will 
be discussed in the following Subsection. In other contexts, when eq. (fT6]l holds, the set Q is 
said to be sampling for the space T-Lm [16] . 

3.2 Exact reconstruction of truncated states 

In our case, there is a convenient way to select the sampling points Zk in such a way that the 
resolution operator A is invertible and explicit formulas for their inverses are available. These 
are given by the set of N points (phases) uniformly distributed on a circumference of radius ^/p^. 

Q = {zk = Vpe'"'^/"^, A: = 0, 1, ... - 1}, (22) 

Denote by 5 = {\zk) ,k = 0, 1, . . . , — 1} the CS subsystem associated with the points in Q 
and by 

= Span(|zo), l^i), . . . , \zN-i)) (23) 

the subspace of Ti spanned by S. For finite A^ we have 'H'^ j^H, so that we cannot reconstruct 
exactly every function ijj € T-L from A^ of its samples = (-^fcl^), but we shall proof that for 
ip GT-Lm we can always provide an exact reconstruction formula. 

Definition 3.1. We define the subspace 1-Lm of truncated functions with a finite number of 
particles M < oo as: 

^A/ = Span(|0),|l),...,|M)). (24) 
The subspace T-Lm is clearly a finite (M + l)-dimensional vector subspace of Ti. 
Theorem 3.2. Given a truncated function ip G T-Lm on C, with a finite expansion 

M 

lijj) = ^ am\m), (25) 

•m=0 
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there exists a reconstruction formula il8\) of i/^ 

N-l 

^{z) = J2 ^k{z)^k, (26) 

fc=0 

from N > M of its samples taken at the sampling points (phases) in i22\). through a "sinc- 
type" kernel given by 

^k{z) = -e(^'=^"--")/2 ^ (zz-'r. (27) 

m=0 

Firstly, we shall introduce some notation and prove some previous lemmas. 



Lemma 3.3. For N > M, the frame operator T : Hm ~^ given by = {{zj^\t/;), zj^ S Q} 
[remember the construction after Eq. U6\}l is such that the resolution operator A = T*T is 
diagonal, A = diag(Ao, • • • , Am); in the basis [24\ ) ofHM, witl^ 

\m{p) = N^p^,m = 0,...,M. (28) 
ml 

Hence, A is invertible inTiM- Therefore, denoting jz^) = A~^\zk), the dual frame, the expression 

N-l N-l 

= ^ \zk){zk\ = ^ \zk){zk\ (29) 

fc=0 k=0 

provides a resolution of the identity inTijif- 

Proof. Using (jlOp . the matix elements of T can be written as: 

Tkn = {Zk\n) = e-P/2^e-2-^"/^ ^ Ay2j-,„ (30) 
vn! 



where J- denotes the Rectangular Fourier Matrix (RFM) [T7] given by = -^^e *27rfcn/7V^ ^ _ 
0, . . . , — 1, n = 0, . . . , M. Then, the matrix elements of the resolution operator turn out to 
be: 

N-l N-l 
Anm = (?T-|.4|m,) = Tk-nTkm = {^n^m)^^'^ J'kn-^km = ^n^nmi (31) 

fc=0 fc=0 

where we have used the well known orthogonality relation for RFM: 

at\^ t* t - \^ ( 2-Ki{n-m)/N\^ _ j N, H {n- m) mod TV = 
JV2^-A„fc-Afcm - Z^y^e j - 1 0, if (n-m)mod7V/0 

fc=0 fc=0 

= ^\n-m)modN,0-, (32) 

and this equals N5n,m N > M. Since A is diagonal with non-zero diagonal elements A„,, then it 
is invertible and a dual frame and a (left) pseudoinverse for T can be constructed, T^"*" = A^^T*, 
providing, according to eq. p^ . a resolution of the identity. ■ 



*The quantities Am are well defined for m G N U {0} and they will be used in the case of an infinite number of 
particles M — <x too. Note also that f{n;p) — \n(p)/N is the density function of the Poisson distribution. 
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Proof of Theorem 13.21 From the resolution of the identity (|29p . any ■0 S T-Lm can be written as 
1^) = YJk=o ^k\zk)-, and therefore ^'(z) = (zjV') = Y,k=o ^k{z\^k)- Using that \zk) = A~^Zk), 
we derive that 



M MM 

kn 



{Azk) = ^ {z\m){m\A ^\n){n\zk)= ^ {Am)-^nLTkn= ^ {Am)Xn^5nmTkr 
m,n=0 m,n=0 m,n=0 

M M -, A/ 

Tkm/N 



m=0 (71=0 m=0 

M ^ M 

^ ^ ^ g-^J/2^mgP/2^-m/2gi2^WA^ = _^{zkZ^,-zz)/2 ^ (^^-1)™ = Hfc(z), (33) 

m=0 m=0 

where we have used the expressions of {z\m) given by eq. ([TO]) and eq. ([22]) . and that p = Zfcifc. 



Remark 3.4. It is interesting to note that eq. (j26p can be interpreted as a sinc-type reconstruction 
formula, where the role of the sine function are played by the functions Hfc(z), satisfying the 
"orthogonality relations" 

EkizO = j^Yl e^-^^'^-')™/^ = Pik, (34) 

m=0 

with P = TTJ^ an orthogonal projector onto a M-dimensional subspace of C''^, the range of T. 
In the case of critical sampling, N = M + 1, the result '^k{zi) = 5ik is recovered (corresponding 
to an interpolation formula), but for the strict oversampling case, > M + 1, a projector is 
obtained to account for the fact that an arbitrary set of overcomplete data ^fc. A: = 0, . . . , iV — 1, 
can be incompatible with \ip) £ Hm-O 

A reconstruction in terms of the Fourier coefficients can be directly obtained by means of 
the (left) pseudoinverse of the frame operator T: 

Corollary 3.5. (Discrete Fourier Transform) The Fourier coefficients am of the expansion |^/^) = 
X^^^o (^ra\fn) of any ip E 'H]\j can be determined in terms of the data "^k = (zkl'ip) ols 

am= , y e^"^^'"/^^fc, m = 0,...,M. (35) 

Proof. Taking the scalar product with (z^l in the expression ()25p of jV'), we arrive at the 
over-determined system of equations 

M 

'^Tkmam = '^k, Tkm = {zk\m), (36) 
m=0 

which can be solved by left multiplying it by the (left) pseudoinverse of T, 7^^ = (7'*7')-i7'* = 
A^^T* . Using the expressions of A~^ = diag(AQ ^, A^^, . . . , A^/), given in Lemma [331 and the 
matrix elements Tfcru given by the formula (jlOp . we arrive at the desired result. ■ 

Remark 3.6. Note that the Fourier coefficients Om are obtained as a (rectangular) discrete Fourier 
transform of the data ^(-Zfc) up to a normalizing factor 1/y/Xm (that is, a modulation by A~^^'^), 
which can be seen as a filter. The expression ([35]) provides a discretization of ([Hj), admitting 
an extension to a fast Fourier Transform. 
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4 Infinite Number of Particles 



In this section we study the more interesting/reahstic case of states with an indeterminate 
(infinite) number of particles. Among these we can find the physical states commented in the 
previous section, for which many of the results presented here simplify due to the finiteness of 
all number operator moments. 

Before studying the partial reconstruction of states from a finite number of phase samples, 
let us introduce some general definitions. 

4.1 Discrete pseudo-frames and pseudo-inverses 

We shall be mainly interested in cases where there are not enough phase samples to completely 
reconstruct a given function ^p. We shall call this case under sampling. However, a partial 
reconstruction is still possible. In these cases, S does not generate a discrete frame but a 
"pseudo-frame", and the resolution operator A would not be invertible. But we can construct 
another operator from T, B = TT* , acting on i'^ . The matrix elements of B are 



therefore B is the discrete reproducing kernel operator, see Eq. ([7|). If the set S is linearly 
independent, the operator B will be invertible and a (right) pseudoinverse can be constructed 
for T, = T*{TT*)^^, in such a way that TT^^ = Ip. As in the previous case there is 
another operator, = T^T acting on % which is an orthogonal projector onto the subspace 
T-L^ spanned by 5. A pseudo-dual frame can be defined as 



Bki = {zk\zi) 



(37) 




(38) 



zi&Q 



providing a resolution of the projector Ps 



Tr^T = ^ \zk){Zk 



Twr= ^ \zk){zk\ = Ps 



(39) 



Using this, a partial reconstruction (an "alias") ^ of V' is obtained, 




(40) 



from its samples ^k = (•^fclV')) through some "Lagrange- like" interpolating functions 



Lk(z) = {z\zk) 



(41) 



fulfilling Lk{zi) = 5ki- The alias is the orthogonal projection of ijj onto the subspace 'H'^ , 
that is, I'i/') = Ps\'4')- The relative (normalized) distance from the exact ip to the reconstructed 
function tj; is given by the relative error function: 



As mentioned above, we shall denote this case by under sampling^ since there are not enough 
data to fully reconstruct tp. In other contexts, a set Q is said to be interpolating if, for an 
arbitrary set of data there exists a S % such that {zk\'ip) = [IS]- This condition is 
satisfied in this case since Lk{zi) = 5ki- 

The two operators A and B are intertwined by the frame operator T, TA = BT- If T were 
invertible, then both A and B would be invertible and 7^,"*" = 7]^ = T^^. This case would 
correspond to critical sampling, where both operators A and B can be used to fully reconstruct 
the function ip. However, in many cases it is not possible to find a set of points Q such that both 
A and B are invertible, that is, there is no critical sampling, or there are not sets Q which are 
sampling and interpolating at the same time. The most common example is the Bargmann-Fock 
space of analytical functions on C, where one can find rectangular lattices which are sampling 
(and therefore A is invertible), or which are interpolating (and thus B is invertible), but not 
both simultaneously [21 [16]. Examples of critical sampling are given by the space of band 
limited functions on M and the set Z, which is both sampling and interpolating, and the space 
of functions on the Riemann sphere (or rather its stereographic projection onto the complex 
plane) with fixed angular momentum s and the set of A^^'^-roots of unity, with = 2s + 1 [17]. 

It should be noted that in the case in which there is a finite number A^ of sampling points z^, 
the space £^ should be substituted by C^, and the operator B can be identified with its matrix 
once a basis has been chosen. 



4.2 Partial reconstruction of states and aliasing 

In the previous Section we have seen that, using A^ phase samples, we can fully reconstruct 
functions ip G Hm for a number of particles up to M = A^ — 1. When the reconstruction of 
a function ip G Ti from a finite number A^ of samples is required, we cannot use the results 
of the previous section since the resolution operator A is no longer invertibl^. However, the 
overlapping kernel operator B turns out to be invertible, and a partial reconstruction can be 
done following the guidelines of the previous Sec. 14.11 (undersampling) . 

The space 7i is infinite-dimensional, and therefore in the partial reconstruction of an arbitrary 
function {ip) = Yl'^=o '^n\n') a considerable error will be committed unless further assumptions 
on the probability amplitudes, a„, are made. Since \ip) is normalizable, the Fourier coefficients 
a„ should decrease to zero, thus even if IV') ^ T^m, if goes to zero fast enough, it will 
"approximately" belong to Hm if the norm of = Y1'^=m+i ^nln) is small compared to the 
norm of \ip), for an appropriately chosen M. Let us formally state these ideas. 

Definition 4.1. Let us define by 

M 

Pm =Y.\'^)('^\ (43) 

m=0 

the projector onto the truncated space Hm for some M < oo. We shall denote by 

^M+lW = - 7-1-T - I |2 • (44) 



'While the operator T : "H — s- C has the same expression as in the previous section, the operator A is an 
infinite dimensional matrix given by: Amn = ^]+pN^]'+qN^n' ' with m — j + pN and n = j' + qN, that is, it is a 
matrix made oi N x N diagonal blocks. 
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the normalized squared distance [similar to \4-S^ ] from 

oo 

\^) =Y,<^n\n) (45) 

n=0 

to its orthogonal projection 

M 

onto the subspace T-Lm- In other words, eM+i('0) sine of the angle between ip and ipM- 

We hope that the (normahzed) error committed when reconstructing ij: from N of its samples 
will be of the order of e]\f{ip) (here M = N — 1), which will be small as long as the Fourier 

coefficients a„ decay fast enough. More precisely, if |a„| < C/n" for some constant C, a > 1/2 

and n > N, then 



°° °° 1 />oo -1 -1 

n4^^ y^-i^'" - 2a-l(iV-l)2«-i' 



(47) 



which says that e^(-0) = 0{jj^^)E 

For the special case in which l^fj) is a physical state, for instance a coherent state i-e. 
when an = e"'^' have that e^(C) = 1 — r(A^, |Cp)/r(A^), which approximately verifies 

for large N: 

eivlCj - I 1 if |^|2 > ^ • (48) 

This means that, taking M = N — 1 and denoting hy p = the mean number of particles, if 
N > p then \Cm) — 10) whereas if < p then \Cm) — 0. This behavior is similar for physical 
states, and guarantees that the partial reconstruction will be rather accurate. See later on Eq. 
(I62p for a related discussion on this "droplet" behavior. 

In the next theorem we provide a partial reconstruction formula for general functions ip £ T-L 
and a bound for the error committed. 

Theorem 4.2. Given ^ gT-L, there exists a partial reconstruction of ^p, in terms of the alias 

N~l 

4>{z) = J2 Lk{z)^k, (49) 

fc=0 

from N of its samples ^k, taken at the sampling points in 1121^) . up to an error 

with i>q{p,N) = (uNy. P^^ ' '^^^ Lagrange-like interpolating functions ^J^^ adopt the 

following form: 

N-l oo 

Lk{z) = -e(^'=-"'=— ATI J2 A,+,^ (zz, i)^-+«^, (51) 

j=0 q=0 



* This condition could be more formally stated by saying that ip belongs to a certain Sobolev space with 
k<a- 1/2. 
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where 

oo 

Xj = Y,^j+qN,j = 0,...,N-l, (52) 

q=0 

are the eigenvalues of the discrete reproducing kernel operator B = TT* [defined in p'T] ) with 
matrix elements B^i = {zk\zi)] and Xn is given by [28\). but now for n = 0,1,2, ... . 

The proof of this theorem is quite involved and we move it to the Appendix [B] in order to 
make the presentation more dynamic. There we also prove an asymptotic expression ()89p for 
vq{p,N), so that the quadratic error (|5U|) approaches zero when N ^ oo, with the asymptotic 
behavior 

4(p, N) < 2e^(V)^l - e%i^P) + 2e%{i;) + 0{N-^), (53) 

as long as p < po{N) given in (j90p . Thus, the reconstruction oi ip hy ip is exact in this limit. 

The Fourier coefficients of the alias ip can also be computed, providing a discrete version of 
the Fourier Transform (1141). which also admits a fast version. 



Corollary 4.3. (Discrete Fourier Transform) The Fourier coefficients On of the expansion |^5[ ) 
can be approximated by the discrete Fourier transform: 

^1/2 , 



)= e'™''/'^*fc, (54) 



Proof. The Fourier coefficients of the alias (j49p are given by: 

Af-l N-\ 
On = {s, nj-i/i) = {s, n\Ps\tl;) = ^ (s, n\zk)'i>k = ^ {s, n\zi){B~^)ki'^k ■ (55) 

k=0 k,l=0 

where Ps is defined in Appendix B. Using that Tin = {zi\s,n) = A„ Tin, given in (f30l) . and the 
expression for the inverse of B given in (j79p . we obtain the final result once the orthogonality 
relation (I32D for Fourier Matrices is used. ■ 



Remark 4.4. The expression for the Fourier coefficients entails a kind of "periodization" of 
the original a^. Indeed, putting 

oo oo oo 

*fc = {Zk\ll^) = ^(Zfc|m)(m|V') = ^ TkmClm = ^ Kl,^Tkniam 
m=0 m=0 m=0 

and using the orthogonality relations ()32p we can write (I54p as: 

^l/2 iV-l oo 
An -M"^ ^ \l/2 7r „ 



1/2 oo 

which implies 



AnmodAf fc=o m=0 
An' 1/2 

■ X] -^j+gA^^i+gJV, j = mod N, (56) 



A„ = A„^p^a„+pAr =^ a„+pAr = W ^^!^±^ , Vp G N . (57) 
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This is the "complex plane" counterpart of the typical aliasing effect for band-unlimited signals 
on the real line.D 

We could think that, for the case eN{tp) = 0, we should recover the results of Section [321 
but we shall see that this is not the case. Before, a process of truncation and filtering of {ip) in 
(j49|) is necessary to recover the reconstruction formula (|26|) for functions "i/; G "Hm (|25|). Indeed, 
if M = — 1, the truncation operation 

M 

iV^Af) = Pm\4^) = ^n\n) (58) 

m=0 

followed by a rescaling (a filter) of the Fourier coefficients 

M c 

^ RI^Pm) = ^an\n) (59) 

renders the reconstruction formula for ^^^(z) = {z\RPm\'^) to the expression (j26p . For arbitrary 
functions ip G 7i, the new bound for the squared normalized error turns out to be 



where we have followed the same steps as in the proof of Theorem 14.21 in Appendix [Bl used that 
the spectral radius p{R) = \\R\\ = 1 + uo{p,N) and that RPmPsPm = Pm- Contrary to (f50]) . 
the new bound (j60p is proportional to e|^('0). If we, moreover, assume a behavior for a„ as in 
(|47l) . then we have that the error ([60]) is of order 0(l/iV2""^). 

Let us comment on an alternative approach to the reconstruction of functions ip for small 
eM+i(V')) which will turn out to be more convenient in a certain limit. Actually, for eM+iii)) « 
1 we have 

ll^-PM^f = 6l,+i(V)||Vf « (61) 

SO that, the reconstruction formula ([261) for ipM = Pm^P would give a good approximation of ip, 
similarly to the approach followed in [18], Section 4. The problem is that, in general, the original 
data = {zk\ip) for ip and the (unknown) "truncated" data '$M,k = {zklPMli^) for Vm are 
different unless {z^IPm = {zk\,yk = 0, . . . , iV — 1, which is equivalent to {zk\PM\zk) = 1,VA: = 
0, . . . , — 1. The following proposition studies the conditions under which such requirement is 
approximately satisfied. 

Proposition 4.5. For large M, the diagonal matrix elements of Pm in l-L^ have the following 
behavior: 

[^ if p > Pc + 

(i.e., an approximation to the Heaviside (unit step) function), where 

Pc = M + l, ac = VM + 1 (63) 
denote a "critical" average number of particles and the standard deviation, respectively. 
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Proof. Using the expression (j30p we have 



M 



M 



Pm{p) = {zk\PM\zk) = V TkmTZk = e ^Y] —.P 



r?i=0 



m=0 



We can compute 



dp 



M 



(64) 



(65) 



Qm{p) is an Erlang-(M + 1) distribution (gamma distribution for integer M) with average 
value of Pc = pQM{p)dp = M + 1 and variance a1 = J^{p — Pc)'^QM{p)dp = M + 1. 
The Erlang-(M + 1) distribution converges to a Normal distribution J\f{pc,cr1) for large values 
of M. Furthermore, we know that the Dirac delta distribution 5 is the limit of gaussians 
5{x — fi) = limo— j>oAA(/i, (7^). This implies that Af{pc,cr'^) — ^^{^^^) for large values of M, 
which means that Pm{p) approaches a Heaviside function for large M. This is the so-called 
droplet in quantum Hall effect language (see e.g. p!9] for an alternative proof in this context). 



Figure m shows a plot of Pm{p) as a function of p for different values of M. It is clear how 
Pm{p) approaches the step function as M grows. 






Figure 1: Pm{p) as a function of p for three different values of M: 10, 100 and 1000. 

Remark 4.6. The matrix elements of Pm in have a circulant matrix structure. In fact, they 
can be seen as a Fourier transform of the coefficients A^. They have the expression {zk\PM\zi) = 
Ck-i{p), where 

^^^P) = ^ E . (66) 

m=0 

Note that Cj^^i = Cf , therefore the only independents elements are Ci ,1 = 0, . . . , y- 

5 Conclusions and Outlook 

Restricting canonical CS states to the (discrete) circle provides an alternative way of analyzing 
functions in the Fock-Bargmann representation. For truncated Hilbert spaces, this finite CS 
subsystems provide finite tight frames with interesting properties. The lack of completeness in 
the general case still allows for partial reconstructions up to a certain error which goes to zero 
when the number of phase samples N ^ oo. We also prove that an almost-exact treatment is 
also possible when the number of phase samples is greater that the mean number of particles 
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p for physical states. An alternative exact approach would be to consider samples of p too, thus 
leading to "circular-lattice" structures. 

Quantum state reconstruction is important in Quantum Optics and Quantum Computation 
and Information. Quantum states are reconstructed by using measurements on an ensemble 
of identical quantum states. In Quantum Optics, phase measurements heavily rely on beam 
splitters and phase shifters to make simultaneous homodyne measurements of the light field 
using photon counting statistics [20] . We feel that our mathematical results might be of some 
physical usage here. 

Our formalism deals with pure states, but it could be easily extended to mixed states, 
allowing to reconstruct accurately quasiprobability distribution functions from discrete phase 
samples. 
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A Phase operator 

A naive way of defining a phase operator O was given by Dirac [5] through the polar decompo- 
sition of the annihilation operator 

a = e^®VAA. (67) 

The inherent difficulties with this approach (later realized by Dirac himself) where clearly 
pointed out by [6] (see also [7]). In fact, one can prove that: 

= e-*®e^® = l-|0)(0|, (68) 

which means that e*® is not unitary, i.e., O is nor hermitian, though it is practically so in a 
macroscopic system {p » 1). Actually, one can see that 

{z\e'^\z) = {z\aM''^/'^\z) = e'^ (l + ^) - (69) 

in a macroscopic system (p >> 1). 

The difficulty of defining a proper (hermitian) phase operator was circumvented in [HI [9]. 
This involves describing a bosonic field mode in a finite (truncated) but arbitrarily large state 
space Hm C of M + 1 dimensions. The phase states are defined as: 



M 

|0) = lim (M + l)-i/2yg*ne|^^^ 

n=0 

where the limiting procedure is necessary in order to normalize the states and it must be taken 
with care (after expectation values are calculated). Working with finite M, one realizes that the 
phase states \Qk) with 

27rA; , , 

ft = «o + -ij-y (n) 
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(for arbitrary ^o) constitutes an orthonormal basis of T-Lm- The hermitian phase operator is 
simply defined by 

M 

Q^Y.^k\0k){0k\- (72) 

fc=0 

These states are, however, unphysical, in the sense that in the hmit when M — )■ (X) their 
mean number of particles goes to infinity 

B Proof of Theorem 14.21 

In this appendix we prove the main theorem 14.21 of our work. Before tackling the proof of this 
theorem, we shall introduce some notation and prove some auxiliary results. 

Lemma B.l. The pseudo-frame operator T ■ H ^ given by 7~(V') = {{zk\'ip),Zk S Q} 
[remember the construction after Eq. U6\) l is such that the overlapping kernel operator B = TT* 
is an N X N Hermitian positive definite invertible matrix, admitting the eig en- decomposition 
B = TDT* , where D = diag(Ao, . . . , \n-i) is a diagonal matrix with Xj given by l[52\) and T is 
the standard Fourier matrix. 

Proof. Let us see that B is diagonalizable and its eigenvalues are given by the expression 
(|52|) . which actually shows that all are strictly positive and hence B is invertible. This can be 
done by taking advantage of the circulant structure of B (see e.g. Appendix B in [17]). Actually, 
using the expression of the CS overlap ([7]), we have: 

Bki = {zk\zi) = e-^'exp (pe^-^(i->^)/^) ^ Q^,, (73) 



where the circulant structure becomes apparent. The eigenvalues of B are easily computed by 
the formula: 

N-l 

Xk = Dkk = {F*BF)kk = e^'"'"/^C^-ne-*2'^™'=/^. (74) 

n,m=0 

If we expand ([73]) as 

oo q 

C; = e-*'j;^e2-'''/^, (75) 

insert this in (j74p and use the general orthogonality relation for Rectangular Fourier Matrices 
(j32|) . we arrive at (|52|). It is evident that > 0, VA: = 0, 1, . . . , A*" — 1, so that B is invertible. ■ 
Following the general guidelines of Subsec. 14.11 we now introduce the projector Pg: 

Lemma B.2. Under the conditions of the previous Lemma, the set {\zk) = '^j^^{S^^)ik\zi) , k = 
0, . . . , A^ — 1} constitutes a dual pseudo-frame for S, the operator = T^^T is an orthogonal 
projector onto the subspace H.^ , where T^^ = T*B~^ is a (right) pseudoinverse for T , and 

N-l N-l 

\zk){zk\ = ^ \zk){zk\ = Ps (76) 

fc=0 A;=0 
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provides a resolution of the projector , whose matrix elements in the orthonormal base of Fock 
states (0) of % exhibit a structure of diagonal N x N blocks: 

Pmn{p,N) = {m\Ps\n) = (AmA„)^/^A„^odAfV-'»)modAr,o, m,n = 0,1,2,..., (77) 



with A„ given by 

Proof. If we define = T*B^^ it is easy to check tliat TT^ = In is tlie identity in C^. In 
the same way, Ps = Tj^T is a projector since Pj = Tj^TTr^T = T^T = Ps and it is orthogonal 
Pg = {T*B~^T)* = T*B~^T = Ps since B is self-adjoint. The resolution of the projector is 
provided by Eq. (j39p . Its matrix elements can be calculated through: 

N-l 

Pmn{p,N) = UliB-'hTkn. (78) 

k,l=0 

The inverse of B can be obtained through the eigen-decomposition: 

N-1 

{B-')ik = {:FD-'T*)ik = ^ E A-ie^-^'^'^-')/^. (79) 

1 /2 

Inserting this last expression and Tkn = A„ J^kn in ([7H|) and using the general orthogonality 
relation ^ for RFM, we finally arrive at ■ 

The matrix elements (|77p will be useful when computing the error function (jSOp for a general 
function (I45p . At some point, we shall be interested in their asymptotic behavior for large N 
(large number of samples). In order to give an explicit expression of this asymptotic behavior 
of Pmnip, N), it will be useful to define the following functions: 

Mp,N) = = Y: . V ^"""' n = 0,...,N-l. (80) 

An ^ (n + uNy. 

u=l ^ ' 

Note that Vn{p^^) = ■^HN,n{p) — 1, where H^^nip) are the generalized hyperbolic functions 
|21j . In terms of i^nip, the matrix elements (177p adopt the following form: 

, , (Aj+,AfAj+tjv)V^ _ „(^+*)^/2 I /8l^ 

E:LoA,+.iv " y(JTW7(JTW)! 1 + -,(P,A^)' ^''^ 

for m = j + sN and n = j + tA^, with j = 0, . . . , N — 1 and s, t = 0, 1, . . . , and zero otherwise. 
In particular, for n,m < N — 1, the projector ()8ip adopts the simple diagonal form 

i'mn(p, A) = — i— — n, m = 0, . . . , A - 1. (82) 

1 + A) 

Let us state and prove an interesting monotony property of I'nip, A). 
Lemma B.3. The functions [8U\) are strictly decreasing sequences of n for p > 0, that is: 

Un{p,N) < Um{p,N) ^ n > m, n, m = 0, . . . , A - 1. (83) 
Proof. The diference between two consecutive terms of the sequence is: 
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u=l 



(n + 1)! 



(n)! 



(n + l + niV)! (n + uN)] 



n = 0,...,N -1. 



Moreover, it can be easily checked that 



(n + 1)! 



< 



(n)! 



-, VuGN,n = 0,...,iV-l 



{n + l+uN)\ {n + uNy. 

which imphes that i^n+i < t^n, n = 0, . . . , N — 1 ■ 

Now we are in conditions to prove our main theorem in this section. 

Proof of Theorem 14. 2t According to ([381) , the pseudo-dual frame is defined by 

Af-l N-l ^ N-l 



1=0 



N 

1=0 j=0 



Thus, the interpolating functions ()4ip read: 

N-l N-l 



N 

1=0 j=0 



Noting that 



= E(^l")('^l^«) = E r=r — 



(84) 



n=0 n=0 

and using the orthogonality relation for Fourier matrices (j32p . we arrive at (j5ip . 

Now it remains to prove the bound (I50p for the error. Decomposing lip) in terms of \ipN-i) 
PN-i\i^) and = (I — P/v_i)|'i/'), we can write 



£^{p,N)Uf = {^P\^P) - {i;N-i\Ps\i^N-i) - {i^h-i\Ps\i^h-i) 
-2Re{iJN-i\Ps\i^N-i)- 



(85) 



Let us start by bounding the term 

N-l N-l 
{tpN-ilPsli^N-i) = E \a'n\'^{n-\Ps\n) = E I 



1 



> -^-^^m', (86) 



n=0 



n=0 



where we have used the expression ()82p . Lemma lB.31 in bounding ^ \ 



> 



1 



0, . . . , N — 1 and the definition ()44p . Next we shall bound the term 

- 2Re{^lJN-l\Ps\i^N-l) < 2\{i;N-i\Ps\i^N-i)\ < 2UN-i\\\\Psi^h-i\ 

< 21 



+y„{p,Ar) - l+uo(p,N)^ 



Vn 



'^N-i\\\\Ps\\Uii-i\\ = 2^1-e%{^)eNmm'', (87) 



where we have used that Ps is an orthogonal projector, therefore its spectral radius is p{Ps) = 
\\Ps\\ = 1, and the Cauchy-Schwarz inequality. Using the same arguments, we can bound the 
remaining term as follows: 



{i;j,_,\Ps\i^^.i) < {i^h-i\Ps\^pN-i) < Un-iWWPsWUn-iW = ^% 
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Putting together all this information in (j85p . we arrive at the bound (|50p for the error. ■ 

It can be easily seen that limjv^oo i^oip,^) = 0,Vp > 0. As a consequence, the error ([50j) 
goes to zero as — )• oo. To obtain the order of magnitude of this error, we shall firstly give an 
asymptotic behavior of fo(p, N) for large A^. 

Proposition B.4. The quantity vq{p^N) has the following asymptotic behavior (as a function of 
N): 

^o(p,Ar) = ^p^ + 0(Ar-2A^), (89) 



as long as 



/{2N)\\n 4_ 1 , . 1 



Proof. Let us see that the first addend {u = 1) of 

oo ^ 
u=l 

is dominant when p < pq{N). Indeed, the quotient between two consecutive terms of the series 
da]) is 

^ ((n + l)Ar)! {2N)V > 

where we have used that the quotient of factorials is decreasing in n G N. If we impose the terms 
of the series (j9T|) to be monotonically decreasing for any u £ N, i.e. 



then we arrive at (pO|) . Thus, the first addend, u = 1, of (|9T]) is the leading term. Using the 
Stirling formula, the asymptotic behavior of the factorial j2N)\ second term, u = 2, gives 

the announced result ([89]) . ■ 
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